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Abstract. We generalize to non-flat geometries the formalism of Simon et al. (2005) 
to reconstruct the dark energy potential. This formalism makes use of quantities similar 
to the Horizon-flow parameters in inflation, can, in principle, be made non-parametric 
and is general enough to be applied outside the simple, single scalar field quintessence. 
Since presently available and forthcoming data do not allow a non-parametric and 
exact reconstruction of the potential, we consider a general parametric description in 
term of Chebyshev polynomials. We then consider present and future measurements 
of H(z), Baryon Acoustic Oscillations surveys and Supernovae type 1A surveys, and 
investigate their constraints on the dark energy potential. We find that, relaxing the 
flatness assumption increases the errors on the reconstructed dark energy evolution but 
does not open up significant degeneracies, provided that a modest prior on geometry 
is imposed. Direct measurements of H(z), such as those provided by BAO surveys, 
are crucially important to constrain the evolution of the dark energy potential and the 
dark energy equation of state, especially for non-trivial deviations from the standard 
ACDM model. 
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1. Introduction 

Recent observations e.g., [H El El HI El El [7] indicate that the present-day energy density 
of the universe is dominated by a "dark energy" component, responsible for the current 
accelerated expansion. 

The leading dark energy candidates are a cosmological constant or a slowly varying 

rolling scalar field e.g.,[Um[T0l[TTl[T2l[T3l[Tl] ^ OT & rev * ew ' though explanations 

in terms of modifications of the Friedman equation ( P51 021 ESI ESI EDI EH E21 E31 E31 125] 
and references therein, under the name of "modifications of gravity" and e.g., [261 E3 
l28l [29] under the name of inhomogeneous models) are also being widely investigated. 

An extended observational effort is being carried out (e.g., SNLS, SDSS, PanStarrs 
etc.) and ambitious plans for the future are being proposed or planned (e.g., DES, 
PAU, BOSS, WFMOS, SNAP, ADEPT, DUNE, SPACE, SKA, LSST) with the goal of 
shedding light on the nature of dark energy. With few exceptions, current constraints 
on the nature of dark energy measure an integrated value over time of the Hubble 
parameter, H(z), which in turn is an integral of its equation of state parameter (w = pj p, 
with p denoting pressure and p density). While these constraints are very tightly 
centered around the cosmological constant value, with a 15% error, the finding that the 
average value of w is consistent with —1 does not exclude the possibility that w varied 
in time e.g., [30]. An emerging technique in dark energy studies uses observations of the 
so-called baryon acoustic oscillations (BAO) [311 [321 [33] . The BAO yield a measurement 
of the sound horizon at recombination, a standard ruler visible at different epochs in 
the lifetime of the Universe: at the last scattering surface through cosmic microwave 
background observations and at lower redshifts through galaxy clustering. In galaxy 
surveys, the BAO scale can be measured both along and perpendicular to the line sight. 
In particular the line-of-sight measurement offers the unique opportunity to measure 
directly H(z), rather than its integral. 

To improve our understanding of dark energy, it is important not only to ask 
whether this dark energy component is dynamical or constant, but also, to constrain 
possible shapes of the dark energy potential. As different theoretical models are 
characterized by different potentials, a reconstruction of the dark energy potential 
from cosmological observations could help discriminating among different, physically 
motivated, models. 

Different approaches to the reconstruction of the dark energy equation of state or 
potential have been proposed in the literature e.g., [3U [351 ESI EH EH [39]. In this paper 
we build upon and generalize the reconstruction technique proposed in [37J to non-flat 
universes. In fact jlQl [H] showed that there can be a degeneracy between geometry and 
dark energy properties, thus analyses to constrain dark energy parameters should be 
carried out varying jointly the geometry of the Universe. 

We then apply this reconstruction to existing determinations of H(z) from ages 
of passively evolving galaxies [37], to new Supernovae data [5], and we forecast the 
constraints on the dark energy potential achievable with the next generation of BAO 
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and Supernova surveys. We find that relaxing the assumption of flatness increases the 
error in the reconstructed dark energy evolution but does not open up degeneracies, 
provided that a modest prior on geometry is imposed: a Gaussian prior on Q k with 
r.m.s. o~k = 0.03. Measurements of H(z) such as those provided by BAO surveys, are 
crucial to constrain the dark energy evolution and to break degeneracies among dark 
energy parameters, for non-trivial deviations from the simplest ACDM model. 

The rest of the paper is organized as follows: in §2 we present our reconstruction of 
the dark energy potential. In §3 we describe the priors used in addition to the present 
and future data sets we consider. Our results on the dark energy potential reconstruction 
are reported in §4. For completeness, we present a similar reconstruction applied to the 
dark energy equation of state parameter as a function of redshift (in §5). We conclude 
in 56. 



2. How to reconstruct the dark energy potential from observations 

Ref. [37] presented a non-parametric method to reconstruct the redshift evolution of the 
potential and kinetic energy densities of the dark energy field, using quantities similar 
to the Horizon- flow parameters in inflation e.g., [121 H3]. As a fully non-parametric 
reconstruction would require the knowledge of H(z) and of H{z), [37J also presented a 
general parameterization of the dark energy potential as a function of redshift, V(z), 
in term of Chebyshev polynomials and showed how to reconstruct the potential as a 
function of the scalar field V(4>) from V(z). In this section we will follow Ref. [37] 
to directly relate the dark energy potential V((j)) with observable quantities, but we 
generalize their approach to non-flat universes. 

We restrict ourselves to classical configurations <fi = <f)(t), that do not break the 
homogeneity and isotropy of space-time. The energy-momentum tensor of this scalar 
field configuration is that of a perfect fluid, with density p<$, and pressure p$ given by 

P0 = K{<j>) + V(<j>) , p^ = K(<f>) - V{<f>) and K = (1) 

where K denotes the kinetic energy of the field. The Friedmann's equations then read: 

R2 = 1 (Pr + P<t> + Pk) , (2) 
a 1 dH 2 ^ 

where n = (or k = 8irG) and pk = —k^, where k is the curvature. In Eq. 

(J3J) we introduced the compact notation p T and pt for the total energy density and 
pressure. Using both Friedmann's equations to solve for the kinetic energy of the field 
and considering a single matter component we obtain: 

ZH\z) - I (1 + z) = k (V ( z) + l - Pm {z) + ~p k (zfj , (4) 

An exact non-parametric reconstruction of V(z) is possible only if H(z) and H{z) are 
known (see Ref. [27J for details). Unfortunately present and near future prospective 
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data do not allow this level of precision. However if V(z) can be parameterized, Eq. 
can be integrated analytically: 



H 2 (z) = (if 2 -|(p m ,o + Pfc,o)) 



:i+z) 



+ % m {z) + Pk (z)) - 2(1 + zf / V(x) (1 + x)- 7 dx . (5) 
<-> Jo 

Hereafter the subscript denotes the quantity evaluated at z = 0. 

An interesting parameterization of the potential involves the Chebyshev 

polynomials, which form a complete set of orthonormal functions on the interval [—1, 1]. 

They also have the interesting property to be the minimax approximating polynomial, 

that is, the approximating polynomial which has the smallest maximum deviation from 

the true function at any given order. We can thus approximate a generic V(z) as 

N 

V(z)~J2 X nT n (x) (6) 

n=0 

where T n denotes the Chebyshev polynomial of order n and we have normalized the 
redshift interval so that x = 2z/z max — 1; z max is the maximum redshift at which 
observations are available and thus x G [—1,1]. Since |T n (x)| < 1 for all n, for most 
applications, an estimate of the error introduced by this approximation is given by Ajv+i- 
With this parameterization, the relevant integral in Eq. (jHJ) becomes: 

/ V(y)(l + y)- 7 dy = ^y2\ n T n {x){a + bx)' 7 dx 

JO 1 „_ n J-l 



n=0 
N 

Zmax 



n=0 

where a = 1 + z max /2 and b = z max /2. These integrals can be solved analytically for 
any order n and F n are known analytic functions. Substituting in ((Sj) we finally obtain: 

N 



H*(z,X i ) = (l + z)*HZ 



1 - 3z "S^ —F (z) 

± <J^ max 7 n\ J 

n=0 



^^m.O ( 1 — Tt — i vT ) — ^k,0 



[i + zfj (i + zy ' (8) 

which relates observable quantities such as the Hubble parameter and Q m ,o with the 
coefficients of the Chebyshev expansion of the potential. Determining the coefficients Aj 
in this manner allows one to reconstruct V(z) as in ([6]). To obtain V(4>), however, 
we would also need to reconstruct 4>(z). This can also be accomplished from the 
determination of the coefficients Aj, through the kinetic energy of the field. From the 
first Friedmann equation we have: 

K{z) = \{!Tz) 2 {l + Z?H2{X - Z) 

= 3k- 1 H 2 {\i, z) - Pm {z) - Pk {z) - V{\i, z) , (9) 
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which can be integrated to obtain (j)(z) and thus V(Xi, <p) from V(aii, z): 



4>(z) - 0(0) = 




z yggggg; z) - 2p m (z) - 2p fc (z) - 2y(\, z) 
(l + z)H(X l , z) 



(10) 



where the ambiguity in sign comes from the quadratic expression for the kinetic energy. 
Typically, if we think of a scalar field rolling slowly along its potential, the plus sign will 
be the relevant one. 

In what follows we will only consider a three parameters model (i.e. N = 2). Even 
with only three parameters, if the fiducial model is a ACDM, the forecasted constraints 
on deviations from a flat potential are rather weak especially at z > 0.5. However, with 
the formulation presented here, one can pose the question of how many dark energy 
parameters are required by the data. Techniques such as cross validation [S] 05] 06] 
could be used to address the issue. 

In the next sections we describe the different observables we consider to probe the 
dark energy potential. 

3. Datasets and priors 

Most probes of dark energy measure integrated quantities, for example Supernovae 
measure the luminosity distance (Il(z). There are two known techniques to reconstruct 
directly H(z). One is through the measurement of the BAO scale in the radial direction 
through dr = c/H(z)dz. Current surveys do not yet have the sufficient statistical power 
to do so and thus current measurements are angle-averaged. However, forthcoming and 
future surveys promise to deliver H(z) determination with % accuracy. BAO surveys, 
of course, also provide measurement of the angular diameter distance oU(z). 

The other technique relies on the measurement of ages of passively evolving galaxies. 
It has been demonstrated by recent observations that massive (L > 2L*) luminous red 
galaxies have formed more than 95% of their stars at redshifts > 4. Since then, stars in 
these galaxies have been evolving passively. They are therefore, excellent cosmic clocks, 
where the age of their stars can be inferred from the integrated stellar spectrum using 
stellar evolution theory. 

BAO surveys along with large samples of type 1A Supernovae are among the leading 
techniques to constrain dark energy and an extensive experimental effort is being carried 
out. For this reason we consider presently available "cosmic clocks" data, presently 
available Supernovae data and future BAO and Type 1A supernova surveys. 

3.1. Priors 

In what follows we assume a flat ACDM with H = 73.2 Km s _1 Mpc _1 and Q m = 0.24 
fiducial model. In all cases we consider Gaussian priors of o H = 8 Km s~ 1 Mpc~ 1 for H , 
°~w m — 0.01 for fl m h 2 and = 0.03 for Q^. This is motivated by the fact that current 
data already constrain these parameters at this level e.g., [47] [2] 0]. We will assume here 
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that these parameters can be constrained at this level by combination of e.g., Cosmic 
Microwave Background experiments (e.g., Planck [48] ) and local determinations of the 
Hubble parameter [49J. 

3.2. Baryon Acoustic Oscillations 

Dark matter overdensities in the early Universe produce acoustic waves in the photon- 
baryon plasma that propagate with the speed of sound until the recombination era, 
when photons decouple from baryons and free stream. The baryon wave then stops 
propagating leaving an imprint at a characteristic distance from the original dark matter 
overdensity: the sound horizon length. This process thus provides a standard ruler at 
which the correlation function of dark matter (and thus of galaxies) should peak (e.g., 
[3T]). Evidence of this peak has already been reported in galaxy surveys e.g., [32, 133] . 
Measuring this standard ruler at different redshifts would provide a powerful probe of 
the expansion history of the Universe and thus of the dark energy potential. 

Baryon Acoustic Oscillations (BAO) can be measured both along and perpendicular 
to the line sight. An angular measurement of the BAO scale at redshift z would then 
give: 

where cIa(z) = 1/(1 + z) f Q c/H(z)dz and rsAo is the BAO scale. This measurements 
of cIa(z) can then be compared to Eqs. (jHJ) or ( |2T| below in §5.) to derive constraints 
on the coefficients Aj or Wi respectively. Alternatively, if the redshift precision of the 
survey is good enough, the BAO scale could be measured along the line of sight as 

Az = H{z)r BAO , (12) 

thus providing a direct measurement of H(z). 

To forecast the errors with which H(z) and cIa{z) will be recovered we make use of 
the formulas derived in [50] , where a grid of BAO surveys was simulated with different 
survey parameters and the accuracy found for the observables was fitted to the following 
formulas: 

a d (Zi) = Xq^\I ^rfnl(Zi) (13) 



where 




M*)=rf~\ hrfM ( 14 ) 



Vi is the volume of the redshift bin z%. The fitting formula was motivated by the 
assumption that the accuracy achievable in the observables will be proportional to the 
fractional error with which the power spectrum can be recovered 



A P / 2 / 1 



P - V N m V nP 



i + rr? , (16) 
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where N m oc Vi is the number of of Fourier modes contributing to the measurement and 
n the number density of galaxies surveyed. Non-linearities tend to erase the acoustic 
peaks via mode-coupling: the function in (115!) takes into account that, at increasing 
redshift, increasingly small scales are in the linear regime. The fitting parameters where 
calibrated on N-body simulations by [50J and found to be: Xq = 0.0148; Xq = 0.0085; 
V = 2.16/Zi 3 ; z m = 1.4; 7 = 0.5. 

Here we will consider two setups that roughly encompass ground-based ( "ground" ) 
and space-based ("space") perspective BAO surveys. The survey parameters are 
summarized in Table [H In both cases we assume that shot noise is unimportant at 
the scale of interest and that the redshift determination is good enough to measure the 
radial BAO signal. Along the way we will also report the results for the angular-only 
BAO. This case will be relevant to photometric surveys that can achieve photometric 
errors better than ~ 4% [5IJ [52] 



Survey 


Area (dg 2 ) 




Z"max 


bins in z 


ground 


10000 


0.1 


1 


9 


space 


30000 


1 


2 


10 



Table 1. Survey parameters of the two BAO surveys considered. 
3.3. Galaxy Ages 

The Hubble parameter depends on the differential age of the Universe as a function 
of redshift via H{z) = dzjdt{\ + z)^ 1 . The feasibility of measuring H(z) from high- 
resolution, high signal-to-noise spectra of passively evolving galaxies was demonstrated 
in [53l [37J EH- Here we use the H(z) determination obtained by [37] and publicly 
available at [55], from a compilation of data at < z < 1.8 and generalize the analysis 
of [37] to non flat universes. Recent studies [SSI E7J have clearly established that massive 
(> 2.2L* ) luminous red galaxies have formed more than 95% of their stars at redshifts 
higher than 4. These galaxies, therefore, form a very uniform population, whose stars 
are evolving passively after the very first short episode of star active star formation 
[56l [571 [58] . Because the stars evolve passively, these massive LRG are excellent cosmic 
clocks, i.e. they provide a direct measurement of dt/dz; the observational evidence 
discards further star formation activity in these galaxies. Dating of the stellar population 
can be achieved by modeling the integrated light of the stellar population using synthetic 
stellar population models, in a similar way to what is done for open and globular clusters 
in the Milky Way. The dating of the stellar population needs to be done on the integrated 
spectrum because individual stars are not resolved and therefore the requirements on the 
observed spectrum are stringent as one needs a very wide wavelength coverage, spectral 
resolution and very high signal-to-noise. Ref. [57] has shown that the spectra of these 
massive LRG at a redshift ~ 0.15 are extremely similar, with differences of only 0.20 
mmag, which is another evidence of the uniformity of the stellar populations in these 
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galaxies. There have been already examples of accurate dating of the stellar populations 
in LRGs ( [59~1 |60"1 [531 137] ) where it has been shown that galaxy spectra with sufficient 
wavelength coverage (the UV region is crucial), wavelength resolution (about 3 A) and 
enough S/N (at least 20 per resolution element of 3 A) can provide sensible constraints 
on cosmological parameters. More details can be found in [371 16T] . 

3.4- Supernovae 

The intrinsic luminosity of Type IA Supernovae (SN) can be accurately predicted from 
the decay rate of the Supernovae brightness. This provides bright standard candles 
that can be observed up to redshifts z > 1. Measuring SN apparent magnitude, the 
luminosity distance can thus be inferred: 

d L {z) = (1 + Z) £ ^rdz . (17) 

If their redshift is also measured, Type IA Supernova provide information on the integral 
of 1/H(z) and hence on the cosmological parameters. In this way Supernovae provided 
the first direct evidence for the accelerated expansion of the Universe [621 [63] . 

We consider present and forecasted Type IA Supernovae data in the analysis of 
Section [U For the present Supernovae data we use the sample of [M]. For future 
Supernovae data we assume 1000 Supernovae distributed in 5 redshift bins between 0.8 
and 1.3 plus a sample of 500 Supernovae at low redshift [651 125]. Table [2] summarizes 
the distribution of the Supernovae considered. 



Mean z 


0.1 


0.85 


0.95 


1.05 


1.15 


1.25 


SN 


500 


231 


219 


200 


183 


167 



Table 2. Redshift distribution of the forecasted Supernovae sample. 

We consider a statistical error on /i = 5\ogdL + K of a^ sta t. = 0.1 due to the 
uncertainty of the corrected apparent magnitudes. We also consider a systematic error 
given by [66] a u , syst . = 0.02(1 + z)/2.7. 

For both Supernovae samples, present and forecasted, we marginalize over the 
absolute magnitude of the sample. 

4. Results 

For the different data sets we compute (or forecast, for future data) the constraints 
on the first three coefficients of the Chebyshev expansion, assuming a flat ACDM with 
Hq = 73.2 Km s~ 1 Mpc _1 and Vt m = 0.24 fiducial model with the priors described above 
and errors on H(z) and oIa(z) and oIl(z) as outlined in §3. For Supernovae and galaxy 
ages data the analysis is performed exploring the likelihood surface via a Markov Chain 
Monte Carlo. For the BAO surveys we use both a Markov Chain Monte Carlo and a 
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Figure 1. 1 and 2 a constraints on the reconstructed potential as a function of redshift 
V(z) from present SN data (top left), SN data from a future- space-based experiment 
(top right), galaxy ages (middle left), "ground" BAO survey (middle right), "space" 
BAO survey (bottom left) and the combination of the two BAO surveys (bottom right). 



Fisher matrix approach, finding good agreement between the two techniques. Here we 
present the results of the Fisher matrix analysis. 

The constraints on the Aj thus derived can be translated into constraints on the 
potential using Eq. ([6]). In Fig. [I] we show the results of this reconstruction for 
Supernovae (present and future, top panels) galaxy ages (left middle panel), "ground" 
BAO survey (middle right), "space" BAO survey (lower left) and the combination of 
the two BAO surveys (lower right). 

Notice that the constraints set by all the datasets considered are strongest between 
z ~ 0.1 — 0.3. This is a consequence of the recent dominance of dark energy in the 
cosmo logical history and translates in a strong linear degeneracy between Ao and Ai as 
discussed below. 
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Figure 2. 1, 2 and 3 a contours for Ao, Ai (upper panels) and Ao, A2 (lower panels) 
from measurements of cLa(,z) (left) and H(z) (right) alone, for a "ground" BAO survey. 
While the fractional error in dA{z) is approximately half of the error in H(z), the direct 
dependence of H(z) on Ai coefficients through Eq. yields stronger constraints on 
these parameters. 



4-1. Interpretation of the reconstructed V(z) 

It can bee seen from Eqs. (1131) and (1141) that the fractional error in d A (z) is approximately 
half of the error in H(z). However, there is a direct dependence of H{z) on \ coefficients 
through Eq. (jSJ), while the relation with dA^z) involves an integral and the bounds 
derived from the information on d A (z) are generally weaker than those derived from 
the measurement of H(z). As an example, in Figs. [2] and [31 we show the 1, 2 and 3 
a constraints that a "ground" or "space" BAO experiment could respectively place on 
Ao, Ai and A2 using only the information on d A (z) (left) and on H(z) (right). Notice 
that the constraints derived from the information on H(z) are much tighter. Indeed, we 
found that, with information on d A (z) or d^z) alone, there is a degeneracy between A 
and A2, as shown in the bottom left panel of Figs. [2] and [3] and the right panel of Fig. 4 
in Ref. [37J. This degeneracy is lifted by data constraining H(z) as can be seen in the 
bottom right panel of Figs. [2] and [3] and the right panel of Figs. 3 and 6 in Ref. [37J . 
This favors spectroscopic surveys, which can measure H(z), over photometric surveys 
with large photo-z errors which can only measure d A (z). 
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Pcrit Pcrit 

Figure 3. 1, 2 and 3 a contours for Xo, Xi and A2 from measurements of oIa(z) 
(left) and H(z) (right) alone, for a "space" BAO survey. While the fractional error in 
oIa{z) is approximately half of the error in H(z), the direct dependence of H(z) on Xi 
coefficients through Eq. yields stronger constraints on these parameters. 



Let us consider more closely the strong degeneracy between A and Ai in the top 
panels of Figs. [2] and [3] and in Figs. 3, 4 and 6 of Ref. [37]. This degeneracy is present in 
all the datasets we considered but it is more pronounced when no information on H(z) 
is available and the sensitivity to the A, coefficients relies in integrals like d^{z) or di{z) 
as for Supernovae data. This degeneracy is described by a linear relation between Ao 
and Ai of the form: 

A = aA! + /5. (18) 

This implies, to first order in the Chebyshev expansion of Eq. (J6]): 

2z 

V(z) = (a- l)X 1 + f3 + X 1 . (19) 

Zmax 

Then, for any value of Ai along this degeneracy, there is a redshift z = z max (l — a)/2 
for which the value of the potential is fixed to V — f3. For all the datasets we found linear 
degeneracies between A and Ai with a < 1 and (3 ~ f^; this means that for all datasets 
the potential is better constrained at low redshift, in order to have V ~ ^a- This reflects 
the fact that cosmological data are more sensitive to the dark energy properties for small 
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z, since at larger redshifts the matter component dominates and the dependence on V(z) 
is subdominant. However the exact transition redshift (where Q m = Qde) depends on 
the shape of the dark energy potential. 

As for the effect of considering non-flat geometries we find that, with the prior 
of o"fc = 0.03 in Qk we consider, the effect on the extraction of the DE properties is 
very small, only slightly increasing the error in the reconstructed parameters. However, 
loosening the prior on can severely spoil the constraints shown here significantly 
worsening the degeneracies among the Aj coefficients. As an example we show in Fig. H] 
the effect on the constraints of the first three Aj by the "space" BAO survey for different 
values of the prior on flfc, 0, 0.03, 0.1 and 1. 

Notice that, in spite of the lack of data for z < 1, the "space" -type survey is placing 
strong constraints also in that redshift region. This is a consequence of a) the priors 
imposed at z = 0, b) the very accurate data for z > 1 and c) the information on w(z) 
enclosed in the (1a(z) constraints, making it possible to interpolate the potential given 
by our smooth parameterization to the low redshift regime. Therefore, the combination 
of the two surveys does not provide a significant improvement of the constraints on V(z) 
over the "space" survey alone, however these constraints are much more robust since 
both experiments now cover the whole redshift range shown in Fig. [TJ 

4-2. Reconstructed V{4>) 

Eq. fllOp also enables us to reconstruct A(f>(z) from Viz) thus constraining V{4>). Fig. [5] 
shows the results of such a reconstruction for the 68% best models for the different 
datasets. Note that, upon integration of Eq. ( flOi) up to z max , a range of A0(z) can be 
obtained up to a maximum value when z = z max . This maximum value will strongly 
depend on the actual model that is integrated and on how strongly the field evolves in 
that model. Thus, not all values for A<j) are allowed and showing the 1 and 2 a contours 
would not be fully correct. Indeed, for the ACDM Eq. ffTO]) will always yield A(p(z) = 
regardless of z max . If the constraints placed by a given data set on the model are tightly 
centered around the ACDM very small values of A<p(z) will be recovered from such 
models. 

5. Reconstruction of the Equation of State 

It is widespread to parameterize dark energy not by the scalar field potential but by 
its equation of state. As long as the equation of state w is > —1, parameterizing 
the dynamics of dark energy through the evolution of its effective equation of state is 
equivalent to considering the redshift evolution of the dark energy potential. However 
if we want to allow w < — 1 then the scalar field description as presented above fails. 

Thus, considering the evolution of an effective equation of state is more general 
than considering the potential of a scalar field. In this case it is easier to relate w(z) 
with the observables. Expanding the redshift dependence of the equation of state in 




Figure 4. 1, 2 and 3 a contours for Ao, Ai and A2 for a "space" BAO survey for 
different values of the prior on flk, 0, 0.03, 0.1 and 1. 
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Figure 5. Instead of reporting the 1 and 2 a contours, we plot the reconstructed V(cj>) 
for the 68% best models for the different datasets (see text for more details). 



Chebyshev polynomials analogously to the expansion of V(z): 



N 



w(z) ~ u)jT{x{z)) 



i=0 



and substituting in the first Friedmann equation we would have: 
H 2 {u u z) ~ Hi [n mfl (l + zf + fi fci0 (l + zf 



N 



+ (1 - Vl mfl - fifc, )(l + ^) 3 exp -Z max y^ y U n G n (z) 



n=0 



where now 



GAz) 
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Ti(x)(a + bx) l dx . 
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Figure 6. 1 and 2 a constraints on the DE equation of state from present SN data (top 
left), SN data from a future- space-based experiment (top right), galaxy ages (middle 
left), "ground" BAO survey (middle right), "space" BAO survey (bottom left) and the 
combination of the two BAO surveys (bottom right). 



Note that in this parameterization the present-day value of w is given by 

N 

w = ^(-l) i u i (23) 

i=0 

Eqs. ( 12014221) are a generalization of sec 3 of [37] to non flat geometries. 

We study the constraints that the datasets above can put on the DE equation of 
state through (12T1) expanding the dark energy equation state up to second order in 
Chebyshev polynomials. We perform forecasts using MCMC's. As in Sec. HJ H , Q m 
and Qk are parameters which we marginalize over. The same priors quoted before are 
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Figure 7. Equation of state parameter for a model lying along the degeneracy given 
by dataset that can only constrain dA or dj,. For comparison also the ACDM case is 
shown (w = —1 ). 

also assumed here. We show the results in Fig. [6j As for the reconstruction of V(z), 
dark energy properties are best constrained at z < 0.3 and the H(z) determination is 
crucial in constraining the dark energy evolution especially for non-trivial deviations 
from a constant equation of state parameter. 

Analogously to the case of the reconstruction of V(z), since the dependence of H(z) 
on w(z) is through an integral, for quantities that depend on integrals of 1/H(z) such 
as oIa{z) or cIl(z), the information on w(z) is even more diluted. This explains the 
weak constraints found in Fig. [6] for all the datasets except the two BAO surveys. The 
deterioration of the bounds comes mainly through a degeneracy between wq and W\ that 
can span down to wq ~ —100 and w\ ~ —100 if only information on (1a or oIl is available. 
This degeneracy is solved with information on H(z). As an example let us consider a 
model lying within this degeneracy. The parameters for this model are: HO = 70.38 
km/s/Mpc, Q m = 0.28, Q k = 0.025, w = -57.51, w x = -82.63, w 2 = -28.31 and its 
equation of state parameter is shown in Fig. [7] along with the ACDM values w — — 1 for 
comparison. From Eq. ( |23l) . we can see that these parameters will still give w = —3.19 
today, however for z = 2, w = —168. 

In Fig. [8] we show the comparison between H(z) and oIa{z) for this model and for 
the ACDM model. From the figure it is clear that information on (Ia(z) alone does 
not suffice to discriminate between the two, while the differences in H(z) between the 
two models are large. Even with these extreme values of the parameters, this model 
mimics the d^z) behaviour of the ACDM, but has a significantly different H(z) which 
oscillates around the ACDM H(z). Thus, upon integrating H(z) to obtain dA(z), the 
regions where the model is above the ACDM compensate the ones where it is below. 
A measurement of H(z), however, can easily distinguish the two models. We were not 
able to reproduce this behaviour to the same degree with a two-parameter description 
of the dark energy equation of state dynamics. This example highlights an important 
open issue in dark energy studies: constraints on dark energy parameters coming from 
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Figure 8. Comparison of H(z) and cIa{z) for the example model described in the text 
and in Fig. [7] and for the A CDM model. 



measurement of integrated quantities depend crucially on the choice of the dark energy 
parameterization [67] : in the absence of a theoretical motivation for a parameterization 
of dark energy properties, forecasts and constraints become crucially model-dependent. 

We have checked that fixing u>2 = and Qk = 0, with just a two parameter 
description of the dynamics of the dark energy equation of state via uo and u>i, we are 
able to recover similar constraints to those found with alternative descriptions with two 
parameters. Note that the linear parameterization in [351 EE] , corresponds to Ui = 
for i > 1, and in particular w = uj — u)\ and w' = 2ui/z max . Finally the linear 
parameterization in a [69, 70] , w = wo + w a z/(l + z) for \w a \ « wo can be closely 
approximated by LO{ = for % > 2, with the constraint (1231 . Ref. [71] pointed out that a 
simple, 2-parameter fit may introduce biases: the expansion (1201) allows one to include 
more parameters by increasing N as the observational data improve. 

The analysis shown here shows that, increasing the number of parameters seriously 
spoils our ability to constrain w(z) except at small redshifts. As for the constraints on 
the potential, the stronger constraints at low redshifts are related to very pronounced 
degeneracies between wq and w±. This degeneracies are much less important for the 
two BAO surveys, but the rest of the datasets considered can only effectively constrain 
the dark energy equation of state at small values of z. For a constant equation of state 
(uji = for i > 0) the bound derived will roughly correspond to the narrowest allowed 
region at small redshifts. 

6. Conclusions 

We have generalized to non-flat geometries the formalism of [37] to reconstruct the 
dark energy potential. This approach makes use of quantities similar to the horizon 
flow parameters used to reconstruct the inflation potential [121113]. The method can, in 
principle, be made non-parametric, but present and forthcoming data do not allow a fully 
non-parametric reconstruction. We have therefore considered a parametric description 
in term of Chebyshev polynomials which, for all our applications, we have truncated to 
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second order. For completeness we have also considered a reconstruction of the dark 
energy equation of state redshift dependence in terms of Chebyshev polynomials, also 
generalizing to non-flat geometries the results of |37J. 

We have considered present measurements of H(z) from ages of passively evolving 
galaxies (53J [37], future Baryon Acoustic Oscillation (BAO) surveys and present and 
future type IA supernova surveys and investigated their constraints on dark energy 
properties. 

We present present and forecast constraints both on V(z) (Fig. [T]) and, more 
interestingly, on V{4>) (Fig. [5]), in sec. 4. Model building for dark energy which rely 
on simple single-field models and provide physically motivated potentials should satisfy 
the constraints shown in the left top and middle panels of Fig. [5j In the future, the 
expected constraints can be as tight as those shown in the two bottom panels of Fig [5j 
More complicated models (multi fields etc.) should produce a redshift evolution of the 
effective dark energy potential which satisfies the constraints in the left top and middle 
panel of Fig. CD The expected future constrains can be as tight as shown in the bottom 
panels of Fig. [IJ 

We find that relaxing the flatness assumption increases slightly the errors on the 
reconstructed dark energy evolution, but does not generate significant degeneracies, 
provided that a modest prior on geometry is imposed = 0.03 (e.g., Fig. HJ). 

Dark energy properties are best constrained at z < 0.3: this is the result of the late- 
time dominance of dark energy. Under the assumptions made here, the most crucial 
being the assumption of a smooth V(z) or w(z), we find that high redshift (z < 2) 
measurements of both H(z) and d& are more powerful than low z measurements. 

When constraining the redshift evolution of both the dark energy potential Viz) 
or the dark energy equation of state parameter w(z) with measurements of integrated 
quantities such as d& or d^, there are large degeneracies among the parameters. These 
degeneracies are greatly reduced or removed with measurements of H(z), such as those 
provided e.g., by future BAO surveys. This is illustrated in Figs. [2] and [3] for the 
potential reconstruction. While the H{z) constraint is generally weaker than the dj\,{z) 
constraints, H{z) is more directly related to the dark energy properties and thus offers 
more powerful dark energy constraints. We have illustrated this with an example of a 
model which lies on the "rfyi-degeneracy" for the reconstruction of w(z). This model 
produces d^(z) and d^z) virtually indistinguishable from that of the ACDM, however 
the H(z) are different and easily distinguishable from BAO measurements with H(z) 
information. 

This highlights an important open issue in dark energy studies: constraints on dark 
energy parameters coming from measurement of integrated quantities such as d& or di 
depend crucially on the choice of the dark energy parameterization [e.g., [67]]: in the 
absence of a theoretical motivation for a parameterization of dark energy properties, 
forecasts and constraints become crucially model-dependent. The dependence of the 
constraints on the assumed dark energy parameterization becomes evident only when 
considering non-trivial deviations from a ACDM model (e.g., deviations from a constant 
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w or generic shape of the potential). This issue is greatly alleviated by measurements 
that carry information on H(z). 
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